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ABSTRACT 

Using a sample of 165 X-ray selected QSOs from seven deep ROSAT fields 
(/(o. 5-2.0 kcV) ~4 x 10~ 15 erg cm~ 2 s _1 ), we investigate the X-ray spectral proper- 
ties of an "average" radio-quiet broad-line QSO as a function of redshift. The QSO 
stacked spectra, in the observers 0.1 — 2 keV band, in five redshift bins over the range 

0. 1 < z < 3.2 apparently harden from an equivalent photon index of T ~ 2.6 at 
z = 0.4 to r ~ 2.1 at z — 2.4 as seen in other QSO samples. In contrast, the spectra 
in the 0.5-2 keV band show no significant variation in spectral index with redshift. 
This suggests the presence of a spectral upturn at low energies (< 0.5 keV). Indeed, 
while at high redshifts (z > 1.0) the single power-law model gives an acceptable fit to 
the data over the full energy band, at lower redshifts the spectra need a second com- 
ponent at low energies, a 'soft excess'. Inclusion of a simple model for the soft excess, 

1. e. a black-body component (kT~ 100 eV), results in a significant improvement to 
the model fit, and yields power-law slopes of T ~ 1.8 — 1.9, for all redshift bins. This 
power-law is not inconsistent, within the error bars, with those of nearby AGN in the 
2—10 keV band, suggesting that the same intrinsic power-law slope may continue 
from 10 keV down to below ~ 0.5 keV. We caution that there is a possibility that 
the spectral upturn observed may not represent a real physical component but could 
be due to co-adding spectra with a large dispersion in spectral indices. Regardless of 
the origin of the soft excess, the average QSO spectrum has important consequences 
for the origin of the X-ray background: the average spectra of the typical, faint, high 
redshift QSO are significantly steeper than the spectrum of the X-ray background 
extending the spectral paradox into the soft 0.1-2 keV X-ray band. 
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1 INTRODUCTION 

Following the discovery of X-ray luminous QSO ~ 25 years 
ago (Lampton et al. 1972), we remain without a detailed 
description of the form and emission mechanism of QSO X- 
ray spectra, despite extensive studies of these objects. 

Satellite missions such as Ginga and EXOSAT indi- 
cated that the X-ray emission of radio-quiet QSO (which 
constitute the vast majority of the QSO population) above 
~ 2 keV can be well described by a power-law with a 
photon index of F ~ 1.9 (Williams et al. 1992, Lawson 
et al. 1992). However, the sources detected by these ex- 
periments which are the X-ray brightest, typically with 
f K > 1 x 10~ 12 erg cm" 2 s"\ nearby (z < 0.2) QSO may 



not be representative of the class as a whole. Reeves et al. 
(1997) observed a small number of bright radio-quiet QSOs 
with ASCA, extending these results to higher redshifts. Re- 
cent long-exposure (~ 100 ksec) ASCA observations (Geor- 
gantopoulos et al. 1997) extend the investigation of the av- 
erage QSO X-ray spectra down to even fainter flux limits of 
~5x 10~ 14 erg cm~ 2 s _1 in the 2-10 keV band. Although, 
the photon statistics are not sufficient to be conclusive, there 
is some evidence that the spectral index of these faint QSOs 
is flat (r = 1.5 ±0.2). 

Observations of the soft X-ray spectrum (< 2 keV) of 
radio quiet, optically selected QSOs made with the Einstein 
observatory show a steeper spectral index (F > 2). This sug- 
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gests that the QSO spectra are concave i.e. they cannot be 
represented by a single power-law over a broad energy band 
(eg Schwartz & Tucker 1988). This steep soft spectral index 
could be interpreted as excess emission (soft excess), at en- 
ergies below ~ 1 keV, above the extrapolation of the X-ray 
power-law slope found at higher energies. Such a soft excess 
is detected in a substantial number of quasars (Masnou et 
al. 1992, Saxton et al. 1993), and in about 50 per cent of 
nearby Seyfert type AGN (Turner & Pounds 1989). Ciliegi 
& Macaccaro (1996) examine the X-ray spectral properties 
of 63 X-ray selected AGN, from the Einstein Medium Sen- 
sitivity Survey, finding T ~ 2.4 with an intrinsic dispertion 
in spectral slope of a ~ 0.4, suggesting that X-ray spec- 
tral slopes are independent of whether selection takes place 
at either X-ray or optical wavelengths. After the launch of 
ROSAT (Triimper 1982), sensitive in the range 0.1 - 2.4 
keV and carrying instruments with improved energy resolu- 
tion over previous missions, information on the soft X-ray 
properties of QSO burgeoned, due, in part, to large surveys 
with flux limits (~ 5 x 10 -15 erg cm _2 s _1 ) orders of mag- 
nitude fainter than those of previous surveys. Schartel et 
al. (1996a) studied the stacked (average) spectra of a sam- 
ple of 908 objects from the LBQS sample of optically se- 
lected QSOs. They find a mean spectral slope of Y ~ 2.5 
and marginal evidence for a flattening of the spectral index 
at higher redshifts. Schartel et al. (1996b) analysed the in- 
dividual spectra of a sample of 55 radio-quiet QSOs from 
the ROSAT all-sky survey. They find again steep spectral 
indices F w 2.5 and no evidence for spectral evolution with 
redshift. Laor al. (1997) also derive a steep spectra index 
r ~ 2.72 ± 0.09 using a sample of 19 radio-quiet, optically 
selected QSOs with z < 0.4. 

The above values for the spectral slope of QSOs consti- 
tute the discrepancy known as the spectral paradox (Boldt 
1987), in that the most common class of X-ray source has 
a spectral index significantly steeper than that of the X- 
ray background (XRB), which has V ~ 1.5 at both hard 
(Marshall et al, 1980, Gendreau et al. 1995) and soft en- 
ergies (Georgantopoulos et al. 1996), and so cannot be the 
major contributor to the XRB flux. A resolution of the spec- 
tral paradox could be the evolution of the mean QSO X-ray 
spectrum with redshift. 

Here, we derive the average QSO spectral index using 
a sample of X-ray selected QSOs detected in seven fields 
from our deep ROSAT survey (Shanks et al. 1991, Geor- 
gantopoulos et al. 1996). Previous work by Almaini et al. 
(1996) in five of our fields finds no variation with redshift in 
the slope of the spectra of QSOs (as determined by analysis 
of hardness ratios in the 0.5 — 2.0 keV band) from five deep 
ROSAT fields. In order to identify the effect of a soft ex- 
cess (emission significant only at £0.5 keV) on the resultant 
power- law slope in the 0.1 — 2.0 keV band, we use here the 
full ROSAT energy band and also perform proper spectral 
fits to the data. In this paper we include two further deep 
ROSAT fields in our investigation, creating a sample of 165 
QSOs covering a broad redshift range, 0.1 < z < 3.2, with a 
relatively high mean redshift (z ~ 1.5). We also include two 
longer expsores (50 ksec each) of the QSF1 and QSF3 fields, 
which, along with the increased range of energies considered, 
results in photon statistics that are effectively doubled over 
those in the Almaini et al. work. Preliminary results from 
this work were presented in Stewart et al. (1994). The goals 
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Table 1. Details of observations used. The observations denoted 
by 1 were made with the PSPC-C instrument, all others were 
made with the PSPC-B. 
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Table 2. Details of redshift bins. 



are to a) detect and quantify the possible evolution of the 
QSO spectral index b) compare the QSO spectra with that 
of the XRB c) examine the possibility for an upturn of the 
average spectral index in the soft band. 



2 X-RAY DATA SELECTION 

We use data taken from deep observations of seven fields 
with the ROSAT PSPC (a summary of the observations 
taken is shown in table ^). Over 300 sources were detected 
within the inner-ring of the detector, at the 4<r level, down 
to a flux of ~ 4 x 10~ 15 erg cm-V^O.S - 2.0 keV). After 
an optical identification program with the AAT, 165 QSOs 
are identified, see Georgantopoulos et al. (1996) for the full 
details of the identification procedure and Shanks et al. (in 
preparation) for the catalogue of sources. Redshifts are in 
the range 0.1 < z < 3.2. Fig. |l| shows the distribution. 

Datasets used in this work are available from the 
LEDAS ROSAT public database. We exclude data from pe- 
riods of high particle background i.e. when the Master Veto 
Rate is above 170 count s" 1 (Plucinsky et al. 1993), typ- 
ically ~ 10 per cent of the data is rejected. Correction of 
the spectra for vignetting and the PSF were peformed with 
the Asterix software package. A background spectrum was 
accumulated from the region inside the inner-ring of the de- 
tector, with all sources detected above the 4a level masked 
out and using a circular region with a radius of 90 arcsec. 

The search for spectral evolution of an average QSO 
requires binning the data into sufficient redshift subsam- 
ples to detect trends over the redshift range. However, one 
must strike a balance between the number of redshift bins 
and the resultant number of objects and photons in each 
bin. Clearly the number of photons must be adequate to al- 
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Redshift (z) 

Figure 1. Redshift distribution, N(z), for the QSO sample, with 
a mean redshift z = 1.44. 



Figure 3. Results from fitting a simple power-law model over the 
energy range 0.1-2.0 keV, with an absorbing column fixed at the 
Galactic value. The upper panel shows the best-fit photon indices 
with lcr errors, the lower panel the reduced \ 2 f° r each redshift 
bin. In the upper panel the best-fit T, from a model fit to stacked 
data from QSOs in all redshift bins, is shown with a dashed line, 
with lcr errors given by the dotted lines. 



Figure 2. Results from fitting a simple power-law model over 
the energy range 0.5-2.0 keV, with absorbing column fixed at the 
mean Galactic value. The upper panel shows the best-fit photon 
indices in each redshift bin, with the dashed line and dotted lines 
showing the best-fit photon index and the lcr error found with 
stacked data from QSOs from all redshift bins. The lower panel 
shows the reduced \ 2 f° r each redshift bin. 



low for full spectral fitting, also, a low number of sources 
per bin may result in the spectrum being dominated by 
one particular source, rather then representing the spectrum 
of a mean source. Details of the redshift bins used, includ- 
ing the total number of background subtracted photons in 
each bin, are given in table ^. In addition to these redshift- 
segregated spectra, we created a spectrum including data 
from all 165 QSOs. We note also that the photons were 
distributed evenly over the fields, thus no one observation 
or field will dominate the results. One of the two observa- 
tions of the QSF3 field, and one of the two observations of 
the QSF1 field was made with the PSPC-C instrument, so 
with these data we use the pspcc_gainl_256.rsp matrix. 
We use the pspcb_gain2_256.rsp matrix for all other ob- 
servations, as they are made after October 1991 with the 
PSPC-B. Spectral fitting is performed for spectra from each 
field and instrument simultaneously. 
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Figure 4. Plot of the co-added spectrum from QSOs, at all 
redshifts, observed with the PSPC-B instrument, with the best- 
fit power-law model (see table ^) shown by the solid line. The 
data/model ratio is given in the lower panel, which exhibits the 
systematic discrepancies typical of our QSO sample when fitted 
with a simple power-law model 



3 SPECTRA AND RESULTS 

Initially we considered the energies in the range 0.5 — 2.0 
keV, and fitted a simple power-law plus absorption column 
model to the data, for comparison with the results from 
Almaini et al. (1996). The Galactic JVh for all our fields 
(Stark et al. 1992, see table jjj) is reasonably low, with a 
range of ~ 1.2 — 3.0 cm -2 and an exposure weighted mean of 
A^h~ 2.0 x 10 20 cm -2 . Hereinafter, unless otherwise stated, 
for the Galactic Nh used in model fitting, we use appropri- 
ate TVh for each of the fields included in each of the spectra. 
The best-fit power-law indices found from the stacked pho- 
tons from all redshift bins is T ~ 2.23 ± 0.07 (table |), with 
a xt value of 1.07 from 157 degrees of freedom, in excel- 
lent agreement with Almaini et al. Errors quoted through- 
out this paper are 90 per cent errors calculated on the basis 
of the conservative assumption that all spectral parameters 
(including normalisation) are interesting. 
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Figure 5. Free-fit absorption, found in the 0.1-2.0 keV band 
using a power-law plus photoelectric absorption model, plotted 
against redshift. We show the results from the PSPC-B spectra, 
and PSPC-C spectra in the uppermost and middle panels respec- 
tively, with la errors. The dashed lines are the mean Galactic val- 
ues appropriate for each of the spectra, from Stark et al. (1992), 
with the dotted lines showing the approximate uncertainty in the 
Stark et al. values. The lower panel gives the reduced x 2 f° r each 
model fit 



Figure 6. Best-fit power-law photon indices (along with la er- 
rors) found over the energies 0.1 — 2.0 keV, when including a red- 
shifted black-body soft-excess component in the spectral model. 
The reduced \ 2 for each redshift bin is given in the lower panel. 
In the upper panel, the dashed line denotes the best-fit photon 
index found when fitting a power-law plus black-body model to 
data from all the QSOs, with la errors shown by the dotted lines. 



Following this we fitted the same model to data from 
each redshift bin individually, resulting in best-fit photon in- 
dices presented graphically in Fig. and in tabular form in 



table |4| A % test shows that this model is formally accept- 
able in each redshift bin, furthermore, we find a non-evolving 
spectral slope, consistent, within errors, in each redshift bin, 
with the value for F found when considering data from all 
redshift bins together. These results are also consistent with 
the Almaini et al. results for the QSO spectral slope as a 
function of redshift [see Fig. 3b of Almaini et al. (1996)]. 

We then tested our data over the range 0.1 — 2.0 keV 
with the same model. The derived photon indices are shown 
in the upper panel of Fig. ^| with the xt given in the lower 
panel. For comparison, the best-fit photon index found for 
data taken from all QSOs is represented by the dashed line 



Figure 7. Best fitting black-body temperatures (top panel) and 
calculated luminosities (lower panel) found for each redshift bin 
whilst fixing the power-law slope to T = 2.0. 



in the upper panel, with the la error shown by the dotted 
lines. Evolution of the derived power-law photon index is 



found, in that F 



(z<1.0) 



> 2.50 and F 



(z>1.0) 



< 2.25. To es- 



tablish the significance of this difference in spectral slope, 
we performed spectral fits on two further spectra, including 
data from QSOs in the redshift range 0.0 — 1.0 and 1.0 — 3.2 
respectively (for results see table H). This difference is sig- 
nificant at over the 8.5a level. It is, however, easily seen that 
this model provides an increasingly inadequate description 
of the data for smaller redshifts, indeed a x 2 test rejects the 
model fit in the two lower redshift bins to at least the 99.9 
per cent level. 

The typical shape of a spectrum that rejects this model 
is shown in Fig. ^, the residuals (lower panel) clearly indicat- 
ing a systematic difference between the data and the model 
at energies J; 0.5 keV. To determine whether the data re- 
quires absorption in excess of the Galactic value, as seen 
in a small number of AGN (Ciliegi & Maccacaro 1996), we 
allowed the absorbing column density to be a free param- 
eter in the fit, again performing the fit for both PSPC-B 
and PSPC-C spectra simultaneously, and with independent 
column densities. Fig. |E] shows how the best fitting column 
density varies with redshift for the PSPC-B spectra (upper- 
most panel) and the PSPC-C spectra (middle panel), with 
the dashed line representing the appropriate weighted mean 
Galactic Ah and the dotted lines the estimated uncertainty 
in the mean. The total absorbing column is consistent with 
the galactic value in all the redshift bins. Neither inclusion 
of Ah as a free parameter nor including an intrinsic warm- 
absorber reduce the \ 2 significantly, as determined by an 
F-test (Bevington & Robinson 1994), in any of the redshift 
bins. 

In an attempt to fit the 0.1 — 2.0 keV spectra, and guided 
by results from nearby AGN which have a spectral upturn 
towards low energies (eg Turner and Pounds, 1989, Fiore 
et al. 1994), we tried a two-component model consisting of 
a power-law plus a redshifted black-body, the derived pho- 
ton indices and xt are shown in Fig. ^| and table ^. We 
see a significant reduction in the \ 2 due to the inclusion 
of this soft-excess even though the model is still formally 
rejected in the lowest redshift bin. However this result is 
probably not unreasonable, as we have co-added the pho- 
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Table 5. Temperature and luminosity of the best-fit black-body 
component. Errors are 90 per cent. The power-law component 
was fixed at T = 2 to allow better constraint of the black-body 
parameters. 



tons from many sources at slightly different redshifts (and 
having different physical parameters and geometries), thus 
a single black-body model may not provide a perfect de- 
scription of the data. We find that the power-law slope de- 
rived when including a black-body soft-excess is, in each 
redshift bin, consistent within the 90 per cent errors with 
r w 1.8 ±0.2, the value found for the entire sample of QSOs 
together. This value is consistent, within the 90 per cent er- 
ror bars, with the intrinsic AGN power-law index of T ~ 1.9 
(Nandra & Pounds 1994), suggesting that the hard power- 
law continues down to below ~ 0.5 keV in the QSO rest 
frame. The best-fit temperatures of the black-body compo- 
nents (kT ~ 100 — 200 eV) are in general agreement with 
those found for quasars (Rachen et al. 1996) and nearby 
radio-quiet QSOs (Fiore et al. 1994). Modelling the soft- 
excess as a thermal bremsstrahlung, a hot plasma, or an 
additional power-law yields similar % 2 results over the five 
redshift bins, and derived power-law indices that, again, are 
around T ~ 1.8. 

While we note that our description of the soft-excesses 
found in the binned spectra as a simple single-temperature 
black-body is simply an expedient to obtain a better statis- 
tical representation of the data and as such may not be fully 
representative of the true spectra of individual objects, or 
the underlying physical properties, it is of interest to exam- 
ine the behaviour of the black-body parameters as a func- 
tion of redshift bin. In Fig. Q and table ^ we present the 
best fit black-body temperature and the luminosity of the 
black-body component as a function of redshift (now fixing 
the power-law photon index at 2 in order to better con- 
strain the black-body properties). It is clear that the black- 
body temperature, while consistent with measurements for 
nearby objects of similar luminosity (Fiore et al. 1994) at 
zero redshift, increases in temperature to higher redshifts, 
rising from ~ 100 eV in the most local bin to of order 200 
eV at redshifts ~ 2. In contrast, the mean luminosity of the 
black-body per object remains approximately constant. 

We first ask whether the effect could be some instru- 
mental artefact, as the black-body temperature in the ob- 
server's frame remains roughly constant. A systematic error 
in the calibration of the effective area, for example, could 
cause such an effect. We can rule out this possibility by not- 
ing that the fractional contribution from the soft excess rel- 
ative to the power-law component varies with redshift from 
bin to bin, while a systematic uncertainty or error in the 
calibration would give a constant fractional contribution. 
Another possibility is that the results we obtain are an arte- 
fact of our assumption that the spectra of QSOs within a 



redshift bin can be represented by a single value of V. For 
example coadding a step and a fiat spctrum would result 
in a steep and flat spectrum at low energies and high en- 
ergies respectively i.e. in a concave spectrum. In principle, 
deriving the spectra of individual QSOs would circumvent 
this problem. Fiore et al. (1994) find that a large fraction 
of the QSOs in their sample (four out of six) reveal a soft 
excess at low PSPC energies. In contrast Laor et al. (1997) 
find that the PSPC spectra of the majority of their LBQS 
QSOs are consistent with a single power-law to within 30%. 
Unfortunately, our objects typically are faint and we can- 
not derive individual spectra. However, as an exercise we 
have derived individually the spectra of the two brightest 
QSOs in the QSF3 field. The column density was fixed at the 
Galactic value N H = 1.7 x 10 20 cm- 2 . For RXJ0342.6-4404 
(z=0.377) we obtain F — 3.18+QQg for a single power- law fit 
(X 2 = 59.9/51 dof). The addition of a black-body compo- 
nent (0.06 keV) is significant at the 99 per cent confidence 
level (x 2 = 51.0/49 dof) with T = 2.58t^- In contrast 
RXJ0342.0-4403 (z=0.635) can be fit with a single power- 
law (r = 2.9+Q io ) without the need of an additional com- 
ponent (x 2 = 40.3/38 dof). It is evident that a conclusive 
answer on the origin of the soft excess cannot be given with 
the photon statistics curently available. Finally, it is pos- 
sible that a range of individual QSO spectra could result 
in artificially induced behaviour which varies with redshift 
due to sampling different energy ranges in each redhsift bin. 
If for example QSOs had a range of spectral indices then 
sampling at different energy ranges would bias the results 
in that using higher energy bands flatter spectra would be 
preferentially be sampled giving rise to an apparent flatten- 
ing of spectrum with redshift. Such an explanation for the 
behaviour observed here seems contrived and unlikely given 
the lack of variation in spectra in bins with redshifts greater 
than ~ 1 . 

Next, we consider the possible physical origin of the 
spectral upturn. The soft excesses observed in AGN spectra 
are often ascribed to the high energy end of the big blue 
bump seen in the optical/UV region of the spectrum and 
thought to originate in the accretion disk (Walter & Fink 
1993, Turner & Pounds 1989). The properties of the accre- 
tion disk are then dominated by the mass accretion rate 
and mass of the central black-hole. In the simplest mod- 
els of the emission from the accretion disk, the spectrum 
is a multi-temperature black-body with the maximum tem- 
perature component coming from the inner disk radii. The 
behaviour of the soft-excess parameters presented here cer- 
tainly suggest that the emission is more complex than would 
be expected in a thin accretion disk. The high temperatures 
found do not fit well with the above models. Moreover, the 
constant luminosity of the soft excess seen in the ROSAT 
band while the power law component increases by a fac- 
tor of 20-30 is more consistent with up-scattering than with 
thermal emission from a disk powered by potential energy 
loss. We note that the constant luminosity of our measured 
soft excesses is similar to that found by Saxton et al. (1993) 
for a sample of nearby QSOs. They found that the soft ex- 
cess luminosity saturated above a power law luminosity of 
~ 10 44 erg s~\ the mean luminosity in our lowest redshift 
bin. They postulated that this was due to a decreasing tem- 
perature of the soft-excess. Such an explanation is contra- 
dicted by our results, which are more in keeping with an 
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Rcdshift Model: PL Model: PL Model: PL+BB 

range 0.5 - 2.0 keV 0.1 - 2.0 keV 0.1 - 2.0 keV 

r xl(y) r x gH r x 2 » 

0.1-3.2 2.23+°;°7 1.07(157) 2.39±°'°2 2.12(231) l-83l°;Jg 1.74(228) 

0.1-1.0 2.32+Jj;^ 1.25(97) 2.60to,ol 2.68(167) 1.80t°'.23 2.09(164) 

1.0-3.2 1.20(108) 2.24±{J;{^ 1.60(168) 1.90±qII 1.54(165) 

Table 3. Details of spectral fitting results. Errors quoted are at the 90 per cent confidence level. 



Redshift Model: PL Model: PL Model: PL+BB 

range 0.5 - 2.0 keV 0.1 - 2.0 keV 0.1 - 2.0 keV 

r xr» r x » r xgM 

to\l 0.65(58) 2.53i Q3 2.23(224) l-84±°;^ 1.34(217) 

+q\ 7 7 1.35(47) 2.54+^;°^ 1.42(159) 1.77t° H 1.12(149) 

-o.il 0.99(60) 2.2llg;°g 1.31(180) 1.57±oll 1.05(169) 

t l 1.23(45) 2.141°;°^ 0.99(128) 2.11±g^| 0.95(117) 

0.89(19) 2.08+°;J° 1.44(47) 2.04+°;^ 1.09(38) 

Table 4. Details of spectral fitting results in five rcdshift bins. Errors quoted are at the 90 per cent confidence level. 
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increasing mass-accretion rate rather than black-hole mass 
being the important factor. 

Regardless of the origin of the spectral upturn at low 
energies, the average QSO spectrum has important implica- 
tions for the origin of the XRB. Boldt et al. (1987) argue 
that as the AGN have hard (2-10 keV) X-ray spectra, signif- 
icantly steeper than the spectrum of the XRB in the same 
band (Marshall et al. 1980), they cannot produce the bulk 
of the XRB. This spectral paradox appears to extend at the 
lower ROSAT energies. The XRB in the 0.5-2 keV band also 
has this flat spectral index, F ~ 1.5, (Georgantopoulos et al. 
1996, Vecchi et al. 1999). All previous studies of the QSO 
spectra in the ROSAT band (eg Schartel et al. 1996, Laor 
et al. 1997) obtain spectral indices much steeper than the 
XRB spectrum in the same band. However, these contain 
mainly bright, low redshift QSOs. In contrast, our sample 
contains typical QSOs in the sense that our objects consti- 
tute a large fraction (over ~ 50 per cent at 1 keV) of the 
XRB and they cover a wide range of redshift and luminosity. 
Our results clearly show that the soft X-ray selected QSOs 
cannot produce the bulk of the soft XRB. Another popu- 
lation with a flat X-ray spectrum is needed at faint fluxes. 
This population may be associated with AGN obscured in 
X-ray wavelenghts; some examples of this population at low 
redshift (z < 0.5) may have been already identified in deep 
ROSAT pointings (Schmidt et al. 1998). 



4 CONCLUSIONS 

We derived the average QSO spectrum and its evolution 
with redshift, using a sample of 165 QSOs from 7 deep 
ROSAT fields (/ x > 4 x 10~ 14 erg cm-V 1 ). We find strong 
evidence (> 8a) for spectral evolution, in that the 0.1 — 2.0 
keV spectra flatten from F ~ 2.6 at low redshifts, to F ~ 2.1 
at z ~ 2.4. At high redshifts a power-law model with Galac- 
tic absorption describes the stacked spectra well, with the 
power-law index in rough agreement with the slope found 
from the 0.5 — 2.0 keV spectra. However, for z< 1.5, the 



spectra of the QSOs in the 0.1 — 2.0 keV band cannot be 
fitted by a single power-law, with Galactic absorption, in 
contrast with some other ROSAT results on low redshift 
QSOs (eg Laor et al. 1997) Furthermore, we find that these 
deviations from a simple power-law cannot be modelled by 
excess Nn absoption along the line of sight. The above sug- 
gest the presence of a soft excess below 0.5 keV in the QSO 
rest-frame. Indeed, when we add a black-body component 
to model the soft excess, we obtain a significant reduction 
in the \ 2 f° r the low redshift bins. However this model is 
still rejected in the lowest redshift bin, probably because a 
single temperature black-body cannot represent adequately 
the coadded soft excess of all QSOs in one bin. The resulting 
power-law component has a photon index F ~ 1.8 ±0.2, and 
shows no evidence of evolution with redshift. This spectral 
index is not inconsistent, within the 90 per cent errors, with 
both the canonical AGN spectral index of F ~ 1.7 observed 
in the 2 — 10 keV band, and the intrinsic spectral slope of 
r ~ 1.9, suggesting that the hard X-ray power-law may con- 
tinue well into the ROSAT band, down to below ~ 0.5 keV. 
The soft excess probably evolves with redshift, maintaining 
a constant luminosity but increasing in temperature from 
100 eV locally, to 200 eV at z ~ 2. We caution that the 
black-body model is simply an expedient to obtain a better 
statistical representation of the data. It is possible that the 
soft excess observed may be just an artefact of the coaddition 
of spectra with a large range of spectral indices. Regardless 
of the origin of the soft excess, the average QSO spectrum 
has important implications for the origin of the XRB. It is 
shown here that the average spectrum of the faint, high red- 
shift QSOs is much steeper than the spectrum of the XRB 
in the soft 0.1-2 keV band, extending the spectral paradox 
to low energies. Future missions such as XMM are expected 
to bring a breakthrough in the spectral studies of high red- 
shift, faint QSOs due to the broad energy coverage and high 
effective area. 



© 0000 RAS, MNRAS 000, 000-000 



The average QSO spectrum and its evolution 7 



5 ACKNOWLEDGMENTS 

AJB acknowledges the receipt of a University of Leices- 
ter studentship. IG and OA acknowledge the support of 
PPARC. 



REFERENCES 

Almaini, O., Shanks, T., Boyle, B. J., Griffiths, R. E., Roche, 
N., Stewart, G. O, Georgantopoulos, I., 1996, MNRAS, 
282, 295 

Bevington, P. R., Robinson, D. K., 1994, Data reduction 
and error analysis for the physical sciences. M c Graw-Hill 
International Editions, New York 

Boldt, E., 1987, Phys. Rep., 146, 215 

Boyle, B. J., Griffiths, R. E., Shanks, T., Stewart, G. C, 

Georgantopoulos, I., 1993, MNRAS, 260, 49 
Ciliegi, P., Maccacaro, T., 1996, MNRAS, 282, 477 
Fiore F., Elvis M., M c Dowell J. C, Siemiginowska A., Wilkes 

B. J., 1994, ApJ, 431, 515 
Gendreau, K., et al., 1995, PASJ, 47, L5 
Georgantopoulos, I., Stewart, G. C, Shanks, T., Griffiths, 

R. E., Boyle, B. J., 1993, MNRAS, 262 619 
Georgantopoulos, I., Stewart, G. C, Shanks, T., Boyle, B. 

J., Griffiths, R. E., 1996, MNRAS, 280, 276 
Georgantopoulos, I., Stewart, G.C., Blair, A.J., Shanks, T., 

Griffiths, R.E., Boyle, B.J., Almaini, O., Roche, N., 1997, 

MNRAS, 297, 203 
Lampton M., Margon, B., Bower, S., Mahoney, W., Ander- 
son, K., 1972, ApJL, 171, 45 
Laor, A., Fiore, F., Elvis, M., Wilkes, B. J., Mcdowell, J. C, 

1997, ApJ, 477, 93 
Lawson, A. J., Turner, M. J. L., Williams, O. R., Stewart, 

G. C, Saxton, R. D., 1992, MNRAS, 259, 743 
Marshall, F. E., et al., 1980, ApJ, 235, 4 
Masnou, J. L., Wilkes, B. J., Elvis, M., Mcdowell, J. C, 

Arnaud, K. A., 1992, A&A, 253, 35 
Matt, G., Fabian, A. C, Ross, R. R., 1993, MNRAS, 264, 

839 

Nandra, K., Pounds, K. A., 1994, MNRAS, 268, 405 
Plucinsky, P. P., Snowden, S. L., Briel, U. G, Hasinger, G, 

Pfeffermann, E., 1993, ApJ, 418, 519 
Rachen J. P., Mannheim K. & Biermann P. L., 1996. A&A, 

310, 371 

Raymond J. C, Smith B. W., 1977, ApJS, 35, 419 

Reeves, J.N., Turner, M.J.L., Ohashi, T., Kii, T., 1997, MN- 
RAS, 292, 468 

Reynolds, C. S., 1997, MNRAS, 286, 513 

Ross, R. R., Fabian, A. G, Mineshige, S., 1992, MNRAS, 
258, 189 

Saxton, R. D., Turner, M. J. L., Williams, O. R., Stewart, 

G. C, Ohashi, T., Kii, T., 1993, MNRAS, 262, 63 
Schartel, N., Green, P.J., Anderson, S.F., Hewett, P.C., 

Foltz, C.B., Margon, B., Brinkmann, W., Fink, H., 

Trumper, J., 1996a, MNRAS, 283 1015 
Schartel, N., Walter, R., Fink, H. H., Truemper, J., 1996b, 

A&A, 307, 33 
Schmidt, M. et al., 1998, A&A, 329, 495 
Schwartz, D.A., & Tucker, W.H., 1988, ApJ, 332, 157 
Shanks T., Georgantopoulos I., Stewart G. G, Pounds K. 

A., Boyle B. J., Griffiths R. E., 1991, Nat, 353, 315 

© 0000 RAS, MNRAS 000, 000-000 



Stark A. A., Gammie C. F., Wilson R. W., Bally J., Linke 
R. A., Heiles C, Hurwitz M., 1992, ApJS, 79, 77 

Trumper, J., 1982, Adv. Space Res., 2, 241 

Turner, T. J., Pounds, K. A., 1989, MNRAS, 240, 833 

Vecchi, A., Molendi, S., Guainazzi. M., Fiore, F., Parmar, 
A.N., 1999, A&A, 349, 73 

Walter, R., Fink, H. H., 1993, A&A, 274, 105 

Wilkes, B. J., Elvis, M., 1987, ApJ, 323, 243 

Williams, O. R., et al., 1992, ApJ, 389, 157 



